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ABSTRACT 



We present updated calculations of stellar evolutionary sequences and detailed 
nucleosynthesis predictions for the brightest asymptotic giant branch (AGB) stars 
/\ '. in the Galaxy with masses between SMq to QM©, with an initial metallicity of 



c3 . Z = .02 ([Fe/H] = 0.14). In our previous studies we used the lVassiliadis fc Wood 



(I1993I ) mass-loss rate, which stays low until the pulsation period r eaches 500 days 



after which point a superwind begins. IVassiliadis fc WoodI ( 119931 ) noted that for 



stars over 2.5Mq the superwind should be delayed until P ^ 750 days at 5Mq. 
We calculate evolutionary sequences where we delay the onset of the superwind 
to pulsation periods of P ~ 700 — 800 days in models of M = 5, 6, and TMq. 
Post-processing nucleosynthesis calculations show that the 6 and 7Mq models 
produce the most Rb, with [Rb/Fe] ~ 1 dex, close to the average of most of the 
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Galactic Rb-rich stars ([Rb/Fe] ^ 1.4 ± 0.8 dex). Changing the rate of the ^^Ne 
+a reactions results in variations of [Rb/Fe] as large as 0.5 dex in models with a 
delayed superwind. The largest enrichment in heavy elements is found for models 
that adopt the NACRE rate of the ^^Ne(a,n)^^Mg reaction. Using this rate allows 
us to best match the composition of most of the Rb-rich stars. A synthetic 
evolution algorithm is then used to remove the remaining envelope resulting in 
final [Rb/Fe] of ~ 1.4 dex although with C/0 ratios > 1. We conclude that 
delaying the superwind may account for the large Rb overabundances observed 
in the brightest metal-rich AGB stars. 

Subject headings: Nuclear Reactions, Nucleosynthesis, Abundances - Stars: Abun- 
dances, Stars: AGB and Post- AGB 



1. Introduction 



Garcia-Hernandez et al.l ( l2006l ) and iGarcia-Hernandez et al.l (120091 ) identified several 



Galactic and Magellanic Cloud intermediate-mass asymptotic giant branch (AGB) stars 
within a sample of OH/IR stars (i.e., bright 0-rich giants with large infrared excesses). 
The large enhancements of the neutron-capture element rubidium (Rb) found by these au- 
thors, combined with the fact that these stars are 0-rich, support the prediction that hot 
bottom burning (HBB) and an efficient third dredge-up (TDU) have occurred in these stars 
( JGarcia-Hernandez et al.ll2007l ). The Rb is hypothesized to be produced by the slow neutron- 
capture process (the s process) in intermediate-mass AGB stars, that in this context have 
masses over M > 4Mq. An enrichment in the element Rb over the elements Sr, Y, and Zr 
is a tantalizing piece of evidence for the efficient operation of the ^^Ne(Q;,n)^^Mg neutron 
source in intermediate- mass AGB stars (JTruran fc Iberu 119771 : ICosner et al.lll980l ). This is 
because the pro duction of Rb is sensitive to the high neutron density ass ociated with this 
neutron source ( JLambert et al.lll995l : lAbia et al.ll200ll : lvan Raai et al.ll2012l ). These observa- 
tional results are particularly important because there is a paucity of observational evidence 
for constraining stellar models of intermediate-mass AGB stars, and especially the efficiency 
of the TDU and HBB. These two last points are much debated in the context of the glob- 
ular cluster abundance anomalies (iKarakas et al. 20021 : IStancliffe et al.l l2004l : iHerwigl |2004J : 
Ventura fc D'Antonall2005at iKarakas et al.ll2006af ). Intermediate- mass AGB stars play an 
important rol e in chemical evolution of galaxies and stella r systems by producing substantial 
nitrogen (e.g.. lRomano et al.ll2010l : iKobayashi et al.ll201ll ). These stars may also prove to be 
imp ortant for the chemica l evolution of a select number of neutron-capture elements such as 
Rb jTravadio et al.ll2004l ). 
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Briefly, during the thermally-pulsing AGB phase the He-burning shell becomes thermally 
unst able every ;^ 10^~^ years for stars with H-exhausted core masses over M^ > 0.8M( 







Herwigl l2005l for a review). The energy provided by the thermal pulse (TP) expands 



(see 

the whole star, pushing the H shell out to cooler regions where it is almost extinguished, 
and subsequently allowing the convective envelope to move inwards (in mass) to regions 
previously mixed by the TP driven convective zone. This inward movement of the convective 
envelope is known as the third dredge- up (TDU), and is responsible for enriching the surface 
in ^^C and other products of He-burning, as well as heavy elements produced by the s 



Busso et al.lll999[ ). In intermediate-mass AGB stars with initial masses > AM^ 



0! 



process (e.g. 

the base of the convective envelope can dip into the top of the H-bu rning shell, causing 
proton- capture nucleosynthesis to occur there (see iLattanzio et al.lll996l for a review). HBB 
nucleosynthesis converts the newly created ^^C into ^^N and can prevent the formation of a 
C-rich atmosphere where C/0 > 1. 



van Raai et al.l ( 12012| ) performed a detailed study on the production of Rb in stellar 



models of intermediate-mass AGB stars covering a range in compositions from solar metal- 
licity down to the metallicity of the Small Magellanic Cloud {Z = 0.004 or [Fe/H0 = -0.7). 
While the qualitative features of the observations could be reproduced (increasing [Rb/Fe] 
ratio with increasing stellar mass or decreasing metallicity) the models could not reproduce 
the stars with the highest Rb abundances, which have [Rb/Fe] > 2.0. A few possible solu- 
tions were suggested including extending the calculations to include models of mass higher 
than 6Mq. Another solution was to consider variations of the mass- loss rate used on the 



AGB phase, which is known to be highly uncertain (e.g., iGroenewegen et al.ll2009l ). 



The objective of our work is to explore two of the possible solutions to the models not 
producing enough Rb: 1) we extend our calculations to include models of up to 9Mq, and 
2) we vary the AGB mass-loss rate in order to increase the number of TPs and TDU mixing 
episodes. For the purposes of the present study we concentrate on models oi Z = 0.02 
(which in this context are slightly super-solar metallicity, see Q. We expect that any of the 
proposed solutions, should they work, would apply equally well to stellar evolution models 
with metallicities appropriate for stars in the Large and Small Magellanic Clouds (hereafter 
LMC and SMC). This is because lower metallicity AGB stars already take many more TPs 
to reach the superwind phase owing to their smaller radii and smaller pulsation periods 
than solar metallicity AGB stars. We also examine the effect of variations in the nuclear 
network and v ariations in the ^^Ne +a reactions, including adopting the latest rates for these 
reactions from llliadis et al.l ( l2010l ). 



^We adopt the standard spectroscopic notation [X/Y] — logio(X/Y)s 
are abundances by number. 



logio(X/Y)Q where X and Y 
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The outline of this paper is as follows. We begin with a description of the new stellar 
evolutionary sequences and in §3] we present the new stellar nucleosynthesis predictions. In 
§l]we present the results of the synthetic evolution to the tip of the AGB and we finish with 
a discussion of the new results obtained in §3 Final remarks are given in §3 



2. New Stellar Evolutionary Sequences 

We calculate the stellar evolut ion and nucleosynthes i s in t wo steps. First, we use the 
stellar evolution code described in iKarakas fc Lattanzid ( 120071 . and references therein) to 
follow the evolution of the stellar structure and abundances important for stellar evolution 
(H, ^He, ^He, ^^C, ^^N, and ^^O) from the zero-aged main sequence (ZAMS) to near the end 
of the AGB phase. Then we perform detailed post-processing nucleosynthesis calculations 
to obtain heavy-element yields from the s tellar evolutionary se quences. In this study, we 
use the same stellar evolution code used in Karakas et al. (2010). which includes the carbon 



and nitrogen-rich low-temperature opacity tables fr omlLederer fc Aringeii (120091 ). The rate 
for the ^'^N(p,7)^^0 reaction has been updated to iBemraerer et al.l (l2()06f). while the rate 
for the triple-a reaction has been updated to NACRE (JAngulo et al.lll999l ). We use the 
mixing length (MLT) theory for convective regions and we set the mixing length parameter 
a = 1.86. The temperature gradient in the convective envelope is determined by the theory 
of convection used in the calculations. It can have a substantial effect on the structure of 
the convective envelope and influences surface properties such as the luminosity and effective 
temperature. It also greatly influences the efficiency of HBB, with the result that a shal- 
lower temper ature gradient results in stro nger burning and higher luminosities (see detailed 
discussion in IVentura fc D'Antonall2005a[). In lower mass AGB s tars, increasing a can also 
increase the efficiency of the TDU (iBoothroyd fc Sackmannlll988l ). In the intermediate- mass 
stars we are considering here, the TDU is already very efficient (A > 0.9, see below) so it 
is not clear that increasing a will lead to a noticeable increase in the production of Rb. 
Regardless, we perform one calculation of a 6Mq, Z = 0.02 model with a = 2.5. 

For this study we concentrate on models with an initial composition oi Z = 0.02, where 
we set the solar metallicity t o be Z(r, = 0.(315, ve ry close to the proto-solar nebula solar 



metallicity of 0.0142 given in lAsplund et al.l (120091 ). The masses considered for this study 
are M = 5, 6, 7, 8 a nd QMc^. For each of these masses we compute one evolutionary sequence 
using the standard IVassiliadis fc WoodI (119931 ) mass-loss prescription. F or the M = 5, 6 an d 



7M(7) m odels we compute one evolutionary sequence using the modified IVassiliadis &: Wood 
( I1993I ) mass-loss prescription described below in §2.11 At this metaUicity, we found that 
M = 8Mq is the limit for producing a C-0 core, and even at this mass we found one 
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carbon fl ash during which a sm all fraction of the C-0 core was burnt to neon (similar to the 
models in lDoherty et al.ll2010l ). The 9Mq model went through off-center carbon ignition and 
(almost) complete core carbon burning and entered the thermally-pulsing AGB with an O- 
Ne core. AGB stars with 0-Ne cores are known as "super- AGB stars" and experience HBB 
with temperatures at the base of the envelope of T > 100 x 10^ K. There are considerable 
uncertainties surrounding the stellar modeling of super- AGB stars including the treatment of 



core c arbon ignition and burning and the efficiency of the TDU (jSiessll2010l : iDoherty et al. 



2010l : IVentura fc D'Antonall201ll ). Here we simply use the 9Mq model as an illustrative 
example of the heavy-element nucleosynthesis in super-AGB stars in comparison to lower- 
mass AGB stars with C-0 cores. 



2.1. Mass loss in bright O-rich AGB stars 

Dealing with the extent and temporal distribution of ma ss-loss in AGB stars is a major 
uncertainty in stellar rnodelin g ( IVentura fc D'Antonall2005bl ). In the existing stellar evolu- 
tion code, the iReimerj (llQTSh mass- l oss pr escription is used on the first giant branch with 
r] = 0.4, and the IVassiliadis fc WoodI (Il993l ) mass-loss prescription is used on the AGB. Ac- 
cording to this prescription the mass-loss rate depends on the radial pulsation period and 
stays low (M < x10~^Mq year~^) until the period, P, exceeds 500 days. After this time 
a luminosity driven superwind begins and the mass-loss rate increases to a few xlO~^Ma 
year~^ and removes the envelope quickly, in a few TPs. 



© 



Vassiliadis fc WoodI (Il993l ) noted that there are optically bright long period variable 
stars with periods of ~ 750 days that are probably intermediate-mass stars of ^ 5Mc7). In or- 



der to prevent their mass-loss prescription from removing these observed objects. IVassiliadis fc Wood 
(Il993l ) recommended a modification to delay the onset of the superwind in stars of masses 
greater than 2.5M0. This suggestion is supported by the optical observations of OH/IR stars, 
which show that the number of obscured (dust enshrouded) stars dranaatical ly increases for 
periods longer than 700 days (see Table 12 in iGarcia- Hernandez et al.l 120071). This modifi- 



cation wBsjT^tjj^ed^ntheAGB model calculations presented by iKarakasI ( I2OIOI ) and used 
in the Ivan Raai et al.l ( 12012| ) study. 



Here we explore the effect of delaying the onset of the su perwind phase on th e AGB 
nucleosynthesis. To do this, we adopt the modification to the IVassiliadis fc WoodI (1993r) 
mass-loss prescription. Specifically, we modify Equation 2 of IVassiliadis fc WoodI ( 119931 ) such 
that the superwind begins at a pulsation period of P ~ 700 — 800 days, depending on the 
initial stellar mass. For comparison to observational data, the IVassiliadis fc WoodI (Il993l ) 



mass-loss prescription gives pre-superwind mass-loss rates between ^ 10 * to ~ 10 ^M( 
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year~^, and superwind mass-loss rates between a few xlO~^M0 year^^ up to a maximum of 
fti 10~^Mq year~^ in the most massive AGB model of 9Mq. These peak values are about a 
factor of three to ten less than the maximum mass-l oss rates measured in red supergiants, 



which have mass-loss rates up to ~ 10 ^Mq year ^ (Ivan Loon et al.lll999al : iDe Beck et al. 



2010l ). iDe Beck et al.l (J2010l ) find a clear relation between the pulsation periods of the AGB 
stars and the derived mass-loss rates, and note that the mass-loss rates level of at 3 x 10~^Mq 
year~^ for periods exceeding about 850 days. This suggests that delaying the onset of the 
superwind to P ~ 700 — 800 days is a suitable although somewhat conservative choice. 



Here in this study we concentrate on variations to the IVassiliadis fc Wood! (119931 ) mass- 
loss presc ription. There are several j ustifications for using the empirical mass- loss prescrip- 
tion from I Vassiliadis &: WoodI (19931) instead of other empirical or theoretical mass-loss pre- 



scrip tions (e.g., 



Arndt et al 



1997 : 



Blocker 1995: Wachter et al. 2002, 2008; van Loon et al. 



20051 ) . Importantly, the model stars we are considering are bright oxygen-rich AGB stars 
(i.e., C/0 < 1) for most of their TP-AGB p hase owing to efficie r it HB B. This means that 
the mass- loss rates of e.g.. lArndt et al.l (119971 ) and lWachter et al.l (120021 ) which were derived 
for C-rich stars, are riot ap propriate for our model stars. In comparison, the sample used by 
Vassiliadis fc WoodI (119931 ) had Galactic and Magellanic Cloud C-rich and 0-rich AGB stars 
and covered a large range of pulsation periods and luminosities. 

Furthermore, the stellar effective temperatures of even the most massive AGB model 
stars we compute are log Teg ~ 3.45. These effect ive temperature s are s imilar to the M-type 
AGB stars in the Hertzsprung- Russell diagram of Ivan Loon et al.l (l2005l . their Fig. 12). This 
means that M-type AGB mass-loss rates are more suitable for our model stars than red 
supergiant mass- loss rates, even though their predicted luminosities are near or just above 
the the cl assical AGB limi t of log L/Lfn ~ 4.75. Using the derived mass-loss rate for M-type 
stars from Ivan Loon et al.l (120051 . in particular, a = —5.65, /3 = 1.05, and 7 = —6.3) and the 



typical stellar effective temperature for the 9Mq star {Tcs ~ 2818K) along with the maximum 
AGB luminosity of L ^ 80, OOOL©, we obtain a mass-loss rate of logM = — 4.IIM0 year^^. 
This is lower than t he ma ximum mass-loss rate of lo gM = — 3.96Mf;^ year~^ f ound using 
Vassiliadis &: WoodI (Il993) for the 9Mq AGB model. iGroenewegen et al.l (120091 ) compared 



the IVassiliadis &: WoodI (119931 ) and Rei mer's mass-loss pres c ription s to a sample of AGB 
mass-loss rates and concluded that the IVassiliadis &: WoodI (119931 ) still provides the best 
description. They did find however that there are also some deficiencies, in particular to the 
maximum adopted mass-loss rates. 
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2.2. Model results 



In Figs. [T] (a)-(c) we show the temporal evolution of (a) the H-exhausted core mass, 
the (b) temperature in the He-shell and (c) the temperature at the base of the convective 
envelope for the SMq, Z = 0.02 models. The model with a delayed superwinci is shown by 



the black solid lines and the model with the standard IVassiliadis &: WoodI ( 119931 ) mass loss is 
shown by the red dashed lines (shifted to the left by t = 1 x 10^ years for clarity). The most 
striking feature of these figures is the extra TPs in the model with a delayed superwind. 
This leads to 1) more TDU episodes and a final higher core mass, 2) more TPs where the 
maximum temperature in the He-intershell exceeds 300 x lO^K and, 3) a longer duration 
for HBB and also a higher peak temperature, which leads to a greater degree of proton- 
capture nucleosynthesis. An examination of Table [1] shows that the features illustrated in 
Fig. [1] are present in all of the stellar models that delay the onset of the superwind. The 
drastic increase in the number of TPs leads to increases in the total amount of He-intershell 
matter mixed into the envelope by the TDU for 5, 6, and 7Mq models. Furthermore, the 
peak temperature in the He-shell increases with TP number and consequently the models 
with a delayed superwind experience more TPs with temperatures sufficient to activate the 
^^Ne neutron source, with important consequences for Rb production as well as other heavy 
elements (see below). 

In Table [T] we show the results of the 6Mq, Z = 0.02 model with a = 2.5 and standard 



Vassiliadis &: WoodI (Il993[ ) mass loss. The strength of the TPs and the efficiency of the 
TDU episodes are almost the same as the model with a = 1.86. Furthermore, the amount 
of He-shell material dredged into the envelope over the AGB phase is almost identical to 
the model with a = 1.86. The most obvious difference is in the efficiency of HBB, which is 
more intense in the model with a higher a. Given the similarities between the two models 
in terms of the amount of He-shell material and the strengths of TPs, we do not consider 
this model further in this study. 

Both the 8Mq and 9Mq models enter the thermally-pulsing AGB with pulsation periods 
P > 500 days (see Fig. |2](a) for the 9Mq model). This means that by the first TP these stars 
are already experiencing the superwind and consequently the behavior of the temperature 
at the base of the convective envelope reaches a peak early on and steadily declines with 
time (see Fig. [2] (d) where the thick solid line represents the temperature during the short 
interpulse periods). The peak in temperature at the base of the envelope is reached in ^ 10 
TPs or so. Similar behavior is noted for the evolution of the bolometric luminosity. For 
example, in Fig. [3] we show the core-mass luminosity relation for the 9Mq model, where the 
luminosity is observed to decrease monotonically with increasing core mass. It is also in- 
structive to examine the behavior of the TDU in these intermediate-mass models. The TDU 



efficiency parameter A is shown to reach a peak at SM© and then decrease with increasing 
stellar mass. This prediction combined with the decreasing He-intershell mass means that 
even though the 9Mq model experienced 163 TPs, the total amount of He-intershell material 
dredged into the envelope is roughly a factor of 2 lower than the 5Mq model with standard 
Vassiliadis fc WoodI (119931 ) mass loss. However, the 9Mq model reaches peak He-shell tem- 
peratures of ~ 400 X lO^K (Fig [2] (c)). This will have a strong impact on the s-process 
nucleosynthesis, and in terms of the final yields, will help to balance out the lower amount 
of dredged-up material. 



Siesd (120101 ) presented the ffist stellar evolution and nucleosynthesis predictions for 
super- AGB s t ars to the end of the AGB phase. Our model compares reasonably well with 
that of ISiesd (120101 ) although the H-exhausted core mass of our model at the ffist TP is 
higher, at 1.17 Mq compared to the same mass and metallicity model computed by Siess (at 
I.O6M0) and that his model experiences 119 TPs compared to our 163 TPs. This could 
be simply a consequence of our model being evolved to a lower final total mass of 3.2M(7Jj, 
at which point Fig. [2] (d) shows that HBB had ceased (the temperature at the base of the 
envelope is ~ 2 x lO^K). Furthermore, our model experiences more efficient HBB with a 
peak temperature of 111 x 10^ K compared to 99.2 x 10^ K in the Siess model, and higher 
temperatures during He-shell burning with a maximum of 402 x 10^ K compared to 341 x 10^ 
K in the Siess model. Finally, of con s equen ce for the production of Rb and other heavy 
elements, the standard models of lSiesd ( I2OIOI ) show no TDU while our 9M(r^ mo del does (see 
Fig. [2] (b) which shows the evolution of the H-exhausted core). ISiessI ( I2OIOI ) calculates a 
series of synthetic post-processing models where he sets the TDU efficiency parameter at a 
constant A = 0.8. We will compare to this model in ^1 . 

Before we move onto discussing the results of the post-processing nucleosynthesis re- 
sults, we briefly me ntion a compari son of the new 5Mq and 6 M(7) models presented he re with 
those presented in iKarakad (120 lOh and used in the previous Ivan Raai et al.l ( 120121 ) study. 
The new stellar evolutionary sequences were computed with updated C and N-rich low- 
temperature opacity tables, which can have the effect of increasing the radius and therefore 
the mass-loss rate when the star becomes C and/or N-rich as a consequence of nucleosyn- 
thesis and mixing (JMarigd l2002() . Furthermore , the newer models use the slower rate for 
the ^^N(p,7)^^0 reaction from iBemmerer et al.l (120061 ) . which effects the rate of H-burning 
during the interpulse phase. For most of the masses considered here, efficient HBB prevents 



^Bascd on information from L. Siess's website, http : //www . astro . ulb . ac . be/^siess/WWWTools/SAGB| 
the final mass of his model was Af « 8.75M0. 

■^Keeping in mind that Siess's synthetic model does not take into account the feedback of the TDU on 
the evolutionary properties of the model star. 
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the formation of a C-rich atmosphere so the question of C enrichment should not be im- 
portant, however all models become N-ric h. At this metallicity, 5Mq is the mass at which 
HBB begins and IVentura fc Marigol (120091 ) showed that the stellar yields of the lowest mass 
models that experience HBB crucially depend on the adopted molecular opacities (see also 
Ventura fc Marigol l2010l ). The new 5Mq model is similar to the previous model (similar 
number of TPs, peak HBB temperatures and luminosities) but dredged up ~ 30% more He- 
intershell material as a consequence of the new model being evolved to a smaller envelope 
mass. The new 6Mq reveals a larger deviation in that it has 15% fewer TP (33TP compared 
to 38TP) and lower peak interpulse luminosities and HBB temperatures by roughly 20% 
and 10%, respectively. This is in spite of the higher mixing- leiigth pa rameter of 1.86 in the 
new models compared to 1.75 previously (IVentura fc D'Antonall2005al . showed that higher a 
leads to higher luminosities and temperatures). The lower temperatures and luminosities is 
likely caused by using the lower ^'^N(p,7)^^0 rate. Overall though, important for s-process 
nucleosynthesis is the amount of material dredged to the stellar surface. Even though the 
new model experienced fewer TPs it was evolved to a lower envelope mass and dredged up 
13% more He-shell material compared to the old model. In summary the newer models, 
while computed with improved input physics, show results that are reasonably consistent 
with the previous calculations. 



3. Nucleosynthesis Results 

The stellar evolutionary sequences described above are fed into a post-processing code 
in order to obtain nucleosynthesis predictions for elements heavie r than iron. The deta ils of 
this procedure have been described in detail elsewhere including iKarakas et al.l (120091 ) and 
Lugaro et al.l (120121 ). Briefly, the nucleosynthesis code needs as input from the stellar evolu- 
tion code variables such as temperature, density, and convective boundaries as a function of 
time and interior mass. The code then re-calculates the abundance changes as a function of 
mass and time using a nuclear netwo rk which cont ains N species and time dependent diffu- 
sive mixing for all convective zones (jCannonlll993l ). For this stud y we utilize thr e e nuc lear 
networks: 1) 166 species using the same rates used in the study of Ivan Raai et al.l (20121 ). 2) 



172 species and based on the latest JINA REACLIB database (ICyburt et al. 



20101 ). and 3) 



a 320 species network for a selection of models to examine predictions for hydrogen through 
to lead. The 166 and 172 species networks include elements from hydrogen to sulfur, and 
then from Fe to Nb (166) or Mo (172) and includes elements of the first s-process peak at 
Rb, Sr, Y, and Zr. 



The production of elements heavier than Fe in intermediate-mass AGB stars crucially 
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dep ends on the rate of t he Ne(a,n) Mg reaction 



the iKarakas et al.l (l2006bl ) rate used by Ivan Raai et al 



I n the 166 species network we use 



mm . whereas in the 172 and 320 



species network we use the NACRE (jAngulo et aLlll999l ) rate. As a further test, we perform 
calculations using the 172 spec ies network and use the latest rates of the ^^Ne(a,n)^^Mg and 



22 



Ne(a,7)^^Mg reactions from llliadis et al.l (120 lOl ). 



F or the initial composition we took the solar distribution of abundances from lAsplund et al 
(120091 ) which sets the global solar metallicity to be Zq = 0.0142 which we round up to 0.015. 
We scale the solar abundances to a metallicity oi Z = 0.02 which results in a slightly super- 
solar initial [Fe/H] of 0.14. Solar abund ances of C, N, O, Ne, Mg, Si, S, Ar, and Fe are the 
pre-solar nebula values from Table 5 of lAsplund et al.l (|2009l ): F is the meteoritic value of 
loge (F)0 = 4.42 from Table 1 of the same paper, chosen because it has a lower uncertainty. 
For the same reason we use the meteoritic values for the solar abundances for many elements 
heavier than Fe including Sr, Eu, and Pb. 

In Figs, m and [5] we show the nucleosynthetic results for the elements lighter and heavier 
than Fe, respectively, from the 5Mq model with a delayed superwind and 320 species. For 
some of the heavier elements, the initial [X/Fe] < (e.g., Sb and Te) because we do not 
include all of the stable isotopes in the network. For example for Te, a large fraction (roughly 
65%) of the elemental component comes from isotopes produced by the rapid neutron capture 
process (i.e., ^^®Te and ^^°Te) and we do not include these isotopes in our 320 species network. 
Overall, the nucleosynthesis is as expected with significant N and Rb production as well as 
some surprises including some Sc production ([Sc/Fe] ~ 0.3 dex), along with C o and Cu 
production (at the level of 0.6 dex for Co and Cu; see also lSmith fc Lambertlll987l ). Carbon 
shows a strong variation with a decrease down to [C/Fe] ~ —0.6 followed by an increase to 
~ 0.07 by the tip of the AGB. Also from Fig. [5] we see no significant production of elements 
beyond the first s-process peak. This is not the case as the met allicity is reduced: t he s- 
process predictions for the 6Mq, Z = 0.0001 model presented by iLugaro et al.l (120121 ) and 
plotted in Fig. [6] shows instead significant Ba and Pb production (final [Ba/Fe] = 1.72 and 
[Pb/Fe] = 0.96). The low-metallicity 6Mq model also produces some Sc ([Sc/Fe] = 0.47) and 
Cu ([Cu/Fe] = 1.0) although little Co. Note that the QMq, Z = 0.0001 model predictions 
illustrated in Fig. would be even more extreme if we were to use the NACRE rate for the 
^^Ne(a,n)^^Mg reaction. 

In Tables [2] and [3] we present a selection of elemental and isotopic abundance ratios at 
the tip of the AGB for the nucleosynthesis models computed for this study. In each case we 
list the mass, the pulsation period (P) at which the superwind begins, the network used for 
the calculation, the rate for the ^^Ne(a,n)^^Mg reaction, and then a selection of isotopic and 
elemental abundance ratios, where we adopt the standard spectroscopic notation, [X/Fe]. 
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All isotopic and element ratios were computed using abundances by number. We find good 
agreement between models using the 172 and 320 species networks and the NACRE rate for 
the ^^Ne neutron source. Table [2] shows that the production of elements beyond the first 
s-process peak (Ba as an example) is negligible and hence our assumption of neglecting these 
species is vahd for this study. Now we compare results from models computed using the 172 



speci es network but different rates for the ^^Ne +a rates (from NACRE and from llliadis et al. 



20101 ) . Using the NACRE rate results in the highest levels of Rb and neutron-capture element 
production although there is only Ri 0.1 — 0.3 dex difference between calculations that use 



the NACRE rate cor npared to those t hat u se the llliadis et al.l (|2010[ ) rate. The 166 species 
network used in the Ivan Raai et al.l ( 120121 ) study results in the lowest levels of neutron- 
capture nucleosynthesis. There are other differences between the two networks besides the 
rate of the ^^Ne +a reactions. The 166 network includes older proton capture rates for the 
CNO cycl e and the NeNa and MgAl chains, along with fits to the neutron-capture cross 
sections of teao et al.l ( 2000 ) whereas the updated 172 network includes JINA REACLIB fits 
to the KADoNiS databasfl 

In terms of the nucleosynthesis of elements lighter than Fe, we find t he effects of HBB are 
clearly evident in all models, including the ^Mq model with a standard IVassiliadis fc Wood 
( 119931 ) mass loss which becomes C-rich at the final TPs. All models show enhanced lithium, 
nitrogen and sodium, mild to moderat e oxygen destruction, and small increases in magnesium 
and aluminum. The observations of iGarcia- Hernandez et al.l ( 120071 ) show that sorne star s 
enhanced in Rb have high Li and some have low Li (see discussion in Ivan Raai et al. l l2012l ). 
The final Li abundances of our models also show a large variation where log e(Li)0 varies 
between 2.4 for the bM^ model with a delayed superwind to —2.2 for the 7Mq with a delayed 
superwind. The 8M0 and 9Mq models have final loge(Li) = 2.0 and 1.9, respectively. 

In all cases the Mg isotopic ratios are also altered from their initial solar values, and most 
of all in the 9Mq model with the most extreme HBB. In this case the alteration of the Mg 
isotopes is mostly the result of proton-captures via the Mg-Al chains. For the lighter elements 
we find good agreement between all nuclear networks for C, N, and O. The [Na/Fe] ratio is 
shown to be significantly lower in the u pdated 172 a nd 32 species networks; this is a result 
of using the newer ^^Na + p rates from lHale et al.l (2004) whereas the 166 species network 



uses the NACRE rates (see discussion in Karakas 



2OIOI ). The elemental [Al/Fe] predictions 



are reasonably consistent between the different networks at a given mass, although the effect 
of using different ^^Ne +a rates is evident in the predictions for the Mg isotopic ratios. 



^Website: http://www.kadonis.orgl 

^Defined by loge(Li) = log(Li/H) +12 and the abundances are by number. 
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In iKarakaa (120101 ) we compared our nucleosynthesis results for AGB stars of < 6Mq 
to other groups. Given the consist ency between the older models and those presented here 
we do not repeat that discussion. ISiesd (120101 ) provides yields - in terms of the mass (in 
Mq) of each element ejected by the stellar winds - for the super- AGB stars considered in his 
study. Here we compare the results from Siess's synthetic 9Mq, Z = 0.02 nucleosynthesis 
model with a constant A = 0.8 and 342 TPs to our 9Mq model, which has A ~ 0.77 and 
163 TPs. For this comparison we use yields from the 172 species network. The hydrogen 
yield is essentially the same in our calculations, our ^He yield is only 3% lower (where the 
final Y Ri 0.38), our ^^C yield is about 60% lower although our yield of ^^N is 3% higher, 
indicating more processing by HBB. We also include a synthetic extension for this model 
(see m]) that results in an extra 60 TPs bringing the total to 223 TPs and increases in ^^C 
but no further increase in ^^N because HBB had ceased. Keeping in mind the qualitative 



nature of the yields presented by ISiessI (|2010[ ) and that our model experienced many fewer 
TPs, these differences are relatively small. 



The study by Ivan Raai et al.l ( 12012| ) investigated the effect of including ^^C-rich regions 
("pockets") in the He intershell of intermediate- mass AGB stars. It is in the ^^C pocket 
that n eutrons are produced via the alternative neutron source, ^^C(a,n)^^0. Ivan Raai et al. 
( I2OI2I ) showed that including ^^C pockets increase the production of Kb by up to 0.4 dex. In 
this study we perform one test calculation where we include some extra mixing of protons 
into the He-intershell of the 5Mq Z = 0.02 model with a delay ed superwind. We d o this 
in the same manner as described by Ivan Raai et al.l (|2012[ ) and iKamath et al.l (|2012[ ) . We 
assume that protons are partially mixed into the intershell over a region of mass, M^ix, at the 
deepest extent of each TDU episode, where the number of protons decreases exponentially 
from the envelope composition, Xp ^ 0.7, to Xp = 1 x 10~^. In this test case we set 
^mix = 1 X 10~^Mq for the 5Mq model, where Mmix is constant. In general, the mass of the 
He-intershell is about a factor of 10 larger than Mmix- The final surface [Rb/Fe] is 0.93 dex, 
an increase of abo u t 0.2 dex on the model without a ^^C pocket. This is smaller than found 
by Ivan Raai et al.l ( 120121 ) at the same ma s s and metall icity, although we no t e tha t the 5Mq, 
Z = 0.02 model used by Ivan Raai et al.l (|2012[ ) with IVassiliadis fc WoodI ( 119931 ) mass loss 
only produced a final [Rb/Fe] = 0.05 dex. These results are encouraging and suggest that 
the inclusion of a small ^'^C pocket could lead to [Rb/Fe] ratios consistent with the average 
of the observed abundances, when considering models with a delayed superwind. 
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4. Synthetic evolution to the tip of the AGB 

In all cases the stellar evolutionary sequences were not evolved until the envelope mass 
was removed and the model star had left the AGB. This is a result of convergence difficulties 
that arise once the envelope mass drops below some value (Mgnv ~ 1^© depending on 
initial mass, see Tabled]). In all cases convergence difficulties end the calculation after the 
cessation of HBB (e.g., see Figs {1] and [2]). The distribution of white dwarf (WD) masses 
reveals the existence of massive WDs (~ 0.8Mq to 1.2Mq) originating from single star 



evolution (JLiebert et al.ll2005t iFerrario et al.l 120051 ). This suggests that whatever the cause 
of the convergence difficulties, real stars somehow lose their envelope at the tip of the AGB. 
The real uncertainty is whether or not the model stars should experience further TPs and 
TDU episodes near the tip of the AGB (JLau et al.ll2012l ). 



Assuming that TPs and TDU episodes will occur, it is possible to estimate the effect 
of these extra TPs on the su rface compositions of our m ode l stars. We d o this in the 



same manner described by e.g.. lKarakas &: Lattanzid (120071) and Ivan Raai et al.l (12012 ) . Our 



method is simple by synthetic AGB evolution standards ( Marigdl200ll : llzzard et al. 



200J) 



alth ough of a similar complexity to the synthetic extensi ons to the n ucleosynthesis used by 



e.g., iForestini fc Charbonnell ( 119971 ) for AGB st ars, and ISiesd (120101) f or sup er- AGB stars. 
For further details we refer to the discussion in iKarakas fc Lattanzid ( 20071 ). Our method 
assumes that all of the physical quantities of the star (e.g., luminosity, interpulse period, 
TE iU efficiency etc) are assumed to remain constant at the tip of the AGB. The study 
by iKarakas &: Lattanzid (120071 ) showed that this is a reasonable assumption at the tip of 
the AGB in intermediate-mass stellar models. The main uncertainty of this method is 
the efficiency of the TDU, which we assume remains the same throughout our synthetic 
calculations. Theoretical models of intermediate-mass stars with deep TD U find that th e 
efficiency does not vary significantly, even with decreasing envelope mass (JKarakasll2010l ). 
For a selection of stellar evolutionary sequences we estimate the number of remaining TPs; 
this is shown in Table |H The stellar mass and pulsation period (P) at which the superwind 
begins are listed together with the amount of material lost at each interpulse period (AM), 
the Rb intershell mass fraction (X(Rb)intersheii), and the final [Rb/Fe] ratio predicted by 
our synthetic calculations. In all cases we take intershell compositions from the 172 species 
network with the NACRE rate for the ^^Ne reaction, as these should provide an upper limit 
to the final [Rb/Fe] ratio. 

For the 9Mq model the synthetic algorithm estimated that there should be 60 TPs 
remaining and the final [Rb/Fe] = 0.90. For this model, the large number of synthetic TPs 
stretches our assumption of a constant luminosity, interpulse period, intershell composition, 
etc. (e.g., see Fig. [3] which shows the evolution of the luminosity with core mass). Varying 
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the interpulse period and luminosity by reasonable amounts (~ 20% over 60 TPs) changes 
the final estimate for [Rb/Fe] by only ^ 0.05 dex. The largest change to the final [Rb/Fe] 
ratio is obtained by increasing the intershell composition of Rb. Over 60 TPs, the intershell 
composition of Rb increased by about a factor of two, as a result of TPs becoming hotter 
with evolution. Assuming that the Rb intershell mass fraction increases again by a factor of 
two we obtain a final [Rb/Fe] of 1.12 dex, which we consider as the maximum value for this 
particular model. Finally, the synthetic evolution algorithm results in a final H-e xhausted 
core mass of 1.19Mq, still well below the critical core mass value of ~ 1.37 Mq (JNomoto 



19841 ). above which electron capture reactions start and the core collapses. This is also the 
case if we were to assume that there are no more TDU episodes, where we predict a final 
core mass of 1.22Mq. 

The main point of this section is to highlight that it is possible to reach [Rb/Fe] values 
of > 1.3 dex in the 6Mq and 7Mq models simply by extending the evolution and allowing 
for TPs near the tip of the AGB. In all of the models, the synthetic extension causes the 
surface composition to become carbon rich, where C/0 > 1, owing to the fact that HBB 
has ceased but the TDU continues. For this reason we conclude that TPs that occur after 
the cessation of HBB are not the cause of the highest Rb overabundances observed in the 
OH/IR stars. 



5. Discussion 



First we examine the behavior of [Rb/Fe] as a function of mass (or bolo metric lu- 
minosity) for the set of models computed using the IVassiliadis fc WoodI (119931 ) mass-loss 
prescription. In Fig. [7] we show the bolometri c luminosity at the t i p of t he AGB versus the 
[Rb/Fe] ratios for models calculated using the IVassiliadis fc WoodI ( 119931 ) mass-loss formula 
and for models with a delayed superwind. We show results from the 172 species network that 
used the NACRE rate for the ^^Ne reaction, and include the final computed and synthetic 
[Rb/Fe] from Tables |3] and HI respectively. The average [Rb/Fe] of 1.4 dex of the Rb-rich stars 
is included along with the maximum uncertainty on the derived abundances. The dashed 
region corresponds to the range of [R b/Fe] where Rb-rich stars are found, where we choose 
a range of 0.6 dex to 2 dex (note that Ivan Raai et al.ll2012l . choose a lower limit of 0.4 dex). 
We can see tha t the most massive AGB stars produce the most Rb, as suggested by the 
observati ons of iGarcia-Hernandez et al.l (20061) and found theoretically for a smaller mass 



range by Ivan Raai et al. 



( 2OI2I ). This trend holds when we compare [Rb/Fe] values obtained 
at the final TP (lower solid line in Fig. [7]) compared to resu lts from the synthetic ex t ension 
(lower dashed line in Fig. [7]). From inspection of Fig. 2 of iGarcia- Hernandez et al.l (|2006[ ) 
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we see that most (12) of Rb-enriched Galactic OH/IR stars have [Rb/Fe] between +0.6 dex 
and +1.6 dex (with an average of +1.4 dex), with a maximum error of ±0.8 dex. Three stars 
have [Rb/Fe] ratios much higher, at [Rb/Fe] ^ 2 ± 0.8 dex. None of the models connected 
by the lower solid line in Fig. [7] produce [Rb/Fe] ratios in the range of the observations, and 
only the synthetic extension can produce high enough [Rb/Fe] ratios in the 8Mq and 9Mq 
models. 

From Fig. [7] we see that models with a delayed superwind produce higher [Rb/Fe] 
abundances (shown by the upper solid line) and are found within the shaded region of the 
observations. This is mostly because these models experience many more TPs and TDU 
episodes. For example, delaying the superwind results in maximum increases of A [Rb/Fe] 
= [Rb/Fe]dsw - [Rb/Fe]vw93 = 0.65, 0.95, and 0.76 dex respectively for the 5, 6, and 7Mq 
models when comparing models that use the NACRE rate. Note that [Rb/Fe] dsw is the final 
surface [Rb/Fe] from the model with a de l ayed s uperwind and [Rb/Fe] vw93 the final [Rb/Fe] 
from the model with IVassiliadis fc Wood! (119931 ). We find smaller increases of A [Rb/Fe] 



0.25, 0. 58, and 0.50 dex when using the 166 species network and the rate from lKarakas et al. 



(J2006bl ). In comparison, we find smaller variations if we compare nucleosynthesis models 
computed with the same stellar evolution input but diff erent nuclear netwo rks and rates 
for the ^^Ne +a reactions. For example, using the newer llliadis et al.l (120 lOf ) rates for the 
^^Ne +a reactions instead of NACRE in the models with a delayed superwind results in 
decreases of up to 0.3 dex in the final [Rb/Fe] abundance. There are larger variations of 
up to 0.5 dex between the 172 species network delayed superwind m odels with the NACRE 
rate and models computed with the 166 species network with the iKarakas et al.l (l2006b[ ) 



rate. Ivan Raai et al.l (120121 ) did not find such large variations and indeed they only show up 
for the models with a delayed superwind. This is a consequence of the models experiencing 
many more TPs at temperatures that allow for activation of the ^^Ne neutron source. 

In Fig. [H] we show the surface [Rb/Fe] abundance versus the pulsation period for two 
models with a delayed superwind. Here the [Rb/Fe] value is seen to increase reasonably 
quickly for pulsation periods higher than about 500 days and 600 days for the 5 and 7Mq 
models, respectively. Once the superwind starts, the [Rb/Fe] abundance is seen to plateau. 
This is because the [Rb/Fe] does not increase very much from one TP to another ([Rb/Fe] 
only increases by ~ 3 — 10 x 10~^ dex per TP). During the superwind the pulsation period 
increases quickly but [Rb/Fe] does not. The predicted [Rb/Fe] of the 7Mq model enters 
the lowest observed abundance range (0.6 dex) of the Rb-rich stars at a pulsation period of 
approximately 700 days, whereas the 5Mq only does so near the end of the calculation. The 
observed period versus [Rb/Fe] relationship for the Rb-rich OH/IR stars is not published 
but a comparison would prove useful for providing extra constraints on the stellar evolution 
calculations. 
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Models with a delayed superwind also produce the highest Rb intershell abundances, 
again a consequence of more TPs. An increase in TPs leads to an increase in core mass, 
which leads to an increase in the maximum He-shell temperature. An increase in temperature 
results in a stronger activation of the ^^Ne(a,n)^^Mg reaction and more Rb production. 
However it is instructive to compare the results for the 9Mq super- AGB model in comparison 
to the 6 and 7Mq models computed with a delayed superwind. These models experienced 
similar numbers of TPs with temperatures over 300 x 10^ K (150, 178, 155 TPs for the 6, 7, 
and 9Mq models respectively) but the 6 and 7Mq models produce mo re Rb in the intershell . 



To understand why we need to examine Rb production in more detail. Ivan Raai et al.l (|2012[ ) 
provided a detailed discussion of the neutron-capture pathways that lead to the production 
of Rb in intermediate-mass AGB stars, here we briefly recap how this happens. Elemental 
Rb consists of two stable isotopes, ^^Rb and ^^Rb. The neutron-rich isotope ^^Rb is produced 
in almost equal quantities in He-intershells of the 6Mq and 9Mq models (at the last TP the 
^^Rb abundance is only about 10% higher in the 6Mq model). However, in the 9Mq model 
the ^^Rb/^^Rb ratio in the He-intershell after the last TP is 0.15 and therefore dominated by 
®^Rb. In comparison the ratio is 0.42 in the 6M(d m odel after the last TP. For comparison. 



the solar ^^Rb/^^ratio = 2.43 ( JAsplund et al.lboOOl ). This indicates that ^^Rb is not being 
produced as efficiently in the 9Mq model. This is because most of the '^^Rb comes from 
the decay of the ground-state of the unstable ^^Kr, which has a half-life of 3939 days (or 
~ 10.7 yearsjj. This lifetime of 10.7 years is longer than the duration of convective TPs 
(r < 3 years) in the 9Mq model, and also less than duration of the TDU (t < 3 years). 
This means that most of the ^^Kr decays to ^^Rb after mixing has finished. Note that the 
next TP is ignited further out in the He-intershell, leaving the previously made ^^Rb in the 
core. In contrast, the duration of convective TPs in the 6Mq model is about 15 years, which 
allows the ^^Kr to decay before the onset of the subsequent TDU episode. 

We can now examine if either of the two proposed solutions we are testing can produce 
enough Rb to match the observations. Increasing the stellar mass does produce more Rb and 
leads to higher [Rb/Fe] ratios. The 8 and 9Mq models can produce enough Rb to match the 
bulk of the Rb-rich stars with [Rb/Fe] ~ 1 dex but only when synthetically extending the 
calculation to account for remaining TPs. Given the uncertainty of the synthetic calculation 
for the 9Mq model we take this result as tentative. Other calculations of s-process nucle- 
osynthesis of AGB stars near the C-0 core limit as well as super-AGB stars are needed to 
verify our synthetic model predictions. The super-AGB models of C. Doherty et al. (2012, 
in preparation) and H. Lau et al. (2012, in preparation) cover a larger range of mass, metal- 



^In our networks we also include as a separate species the meta-stable state of ^^Kr which has a half-life 
of 4.48 hours and also /3-decays to *^Rb. 
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licities, and mass-loss rates than we have explored here. It will be interesting to see if these 
models will be able to produce [Rb/Fe] abundances of > 2 at solar metallicity. 

We can now also address the question as to whether delaying the superwind can produce 
enough Rb to match the [Rb/Fe] ratios derived for the Galactic OH/IR stars. Without 
the synthetic extension, the 6Mq and 7Mq models with a delayed superwind produce the 
highest [Rb/Fe] ratios of 1.14 and 1.02 dex, respectively. Both of these models are within 
the observed range of most Rb-rich stars which have [Rb/Fe] < 1.6 ±0.8 dex. The synthetic 
extension to the nucleosynthesis calculations results in [Rb/Fe] ratios of 1.3 and 1.34 dex for 
the 7Mq and 6Mq models, respectively. This means that the 6Mq with a delayed superwind 
can match the composition of the two Rb-enriched stars that have [Rb/Fe] ~ 2 dex after 
allowing for the substantial uncertainties of ±0.8 dex (only one other star has a higher ratio 
at [Rb/Fe] = 2.5 dex). Overall, we conclude that delaying the superwind leads to higher 
[Rb/Fe] ratios, with most of the Rb-enriched stars (with [Rb/Fe] ^1.6 dex) explained by 
the new stellar evolution models, even without allowing for the synthetic extension. It is 
also satisfying t hat this simple solution i s motivated by observations of intermediate-mass 
AGB stars (e.g. IVassiliadis fc Woodlll993uGarcia-Hernandez et al.ll2007| ). 



The synthetic extension (dashed lines in Fig. [7]) that accounts for final remaining TPs 
after the end of HBB allows us to reach [Rb/Fe] abundances of ~ 1.4 dex. The increase in 
heavy elements comes with a substantial increase in ^^C and in all cases the final predicted 
C/0 > 1. All our model stars would be observ ed as dust-enshrouded bright C-rich AGB 
stars (as seen in the LMC, Ivan Loon et al.lll999bl ). In all cases the synthetic extension takes 
place when the pulsation period is very large (see Fig. [8]) and the strong variability present 
on the AGB has almost ceased. Interestingly, many dust-enshrouded early post-AGB stars 
(with no signs of strong variability or very large periods) are found to show an inner C-rich 
dust shell together with an outer 0-rich dust envelope, whic h is consistent with a late change 
of the chemistry of the central star (see lBunzel et al.ll2009l ). These authors concluded that 
because all of the post-AGB stars still host masers in an 0-rich environment the conversion 
of the inner shell chemistry cannot have happened more than a few thousand years ago. All 
stars must have experienced the conversion into a C star while they were losing the last of 
their envelopes before exposing their cores. 

For the 6 and 7Mq models, only 4 and 7 extra TPs and TDU episodes are needed to 
reach [Rb/Fe] abundances of 1.4 dex. These synthetic estimates are produced with a much 
reduced envelope mass owing to the effect of the luminosity-driven superwind. The smaller 
dilution leads to larger [Rb/Fe] more quickly than during the normal evolution but it also 
perhaps suggests that we are too modest with our delay in the onset of the superwind. For 
some stars the superwind may not begin until the pulsation period reaches P > 1000 days 
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or higher and this could lead to [Rb/Fe] ^1.4 dex. This explanation could also account for 
the composition of the three most Rb-enriched stars, with [Rb/Fe] > 2 dex. 

We studied the inclusion of a ^^C-rich region in the He-intershell of a SM©, Z = 0.02 
model and found it to increase the [Rb/Fe] ratio by ~ 0.2 dex. This suggests that we may 
be able to obtain the average observed Rb value by 1) delaying the onset of the superwind 
even further than done here, and/or 2) including a ^^C pocket. There are questions about 
the formation of ^^C pockets in intermediate-mass AGB models. This is because the tem- 
perature at the base of the conv ective envelope during the TDU ma y become hot enough for 
proton-capture nucleosynthesis (JGoriely fc Siessll2004l : lHerwig||2004l ). For thi s study the most 
impo rtant consequen ce is the l ikely inhibition of formation of the ^^C pocket (JGoriely &: Siess 
20041 ). The models of lHerwid ( 120041 ) have very deep TDU into the CO co re, as a conseq uence 
of the diffusive convective overshoot scheme that he employs (see, e.g.. iHerwigl l2000f ) . One 
interesting consequence of this overshoot is that each subsequent TP will be ignited deeper 
in the core and will be hotter. This would lead to a stronger activation of the ^^Ne(a,n)^^Mg 
reaction, and more Rb production. Also, the very deep dredge-up could reach the Rb left 
in the core, further enhancing the surface composition. It remains to be seen if this can 
produce Rb compositions as high as observed in the OH/IR stars. 



While it has been assume d that the metallicities of the Galactic OH/IR stars are solar 
( IGarcia-Hernandez et al.ll2007 ) , it ca nnot be discounted that there may be a small spread in 
metallicity. In Ivan Raai et al.l (20121) and from Fig. IHwe see that a reduction in the initial 
metallicity results in a substantial increase in the productio n of neutron-cap ture elements 
and of Rb, in particular. For exa mple, using models from iKarakad (120 10[ ) and the 172 
species network that uses the latest llliadis et al.l (120101 ) rate, we find that a GMq, Z = 0.008 



([Fe/H] =-0.3) model produces [Rb/Fe] = 1.05 dex, and a 4Mq, Z = 0.004 ([Fe/H] = -0.7) 
model makes [Rb/Fe] = 1.02 dex. While these [Rb/Fe] abundances are much lower than 
the derived [ Rb/Fe] values of the Rb-rich st ars in the Large and Small Magellanic Cloud 
(up to 5 dex, iGarcia- Hernandez et al.l 120091 ) . they are within the abundances derived for 
the Galactic OH/IR stars. Perhaps the OH/IR stars with the highest [Rb/Fe] abundances 
of > 2 evolved from slightly lower metallicity AGB stars of 6-7Mq, [Fe/H] > —0.3 and a 
delayed superwind? Observational evidence for this hypothesis comes from the metallicity 
distribution of Galactic planetary nebu lae (PNe) that span a rang e from 0.1 — 2.0Zq, ac- 
cording to their O and Ar abundances (ISterling fc Dinersteinll2008l ). Note that this is also 
the case for Type I PNe, which are known to be the descend ants of more massive AGB stars, 
according to their spatial distribution and kinematics (e.g., ICorradi fc Schwarzlll995[ ). Also, 
a lar ge spread of metallicitie s is known in other galaxies such as the Large Magellanic Cloud 
(e.g., iGrocholski et al.ll2006[ ). In summary, while the solutions proposed to explain the three 
most Rb-enriched stars seem promising, further detailed calculations are needed to test these 
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scenarios. 



In Ivan Raai et al.l (|2012[ ) we provided a detailed discussion of the serious uncertainties 
affecting the model atmospheres used to derive abundances for the OH/IR AGB stars, in- 
cluding NLTE and 3D effects on the formation of Rb lines in the atmospheres. Effects caused 
by using ID model atmospheres compared to 3D are expected to be lower in the Galactic 
sources than in the Magellanic Cloud AGB stars owing to their higher metallicities. This 
means that our new results may bring observations and models into agreement without hav- 
ing to invoke large NLTE corrections. Note also that the NLTE corrections are expected to 
be larger in the outliers (the three objects with [Rb/Fe ] > 2 dex), because th e lines are more 



saturated. Finally, we do not repeat the discussion in Ivan Raai et al.l (120121 ) relating to the 
problem of the [Rb/Zr] ratio in the model stars compared to the observations. The problem 
stems from the fact that the models that produce Rb also produce some Zr, whereas Zr is 
absent from the atmosphere of the observed AGB stars. Here we note that the new models 
produce higher [Rb/Zr] ratios, in better agreement with the observations, but it still seems 
likely that an important fraction of the freshly manufactured Zr condenses into dust. 



6. Concluding remarks 

We have shown that delaying the superwind in intermediate-mass AGB stars results in 
far greater production of Rb and other neutron-capture elements at solar metallicity. This 
indicates a possible solution to the discrepancy between Rb observations and theoretical 
model predictions. Including a synthetic model to account for the final remaining TPs 
further increases the Rb composition of the envelope. The extra Rb comes with an increase 
of C and the final C/0 ratio is > 1 for all models. However this exercise shows that not 
many more TPs are needed for the 6Mq and 7Mq models to match the average [Rb/Fe] 
ratio of most of the Rb-rich OH/IR stars. This could be achieved by delaying the start of 
the superwind mass-loss phase even further than done here or including a small ^'^C pocket. 
These solutions may also allow us to reach the [Rb/Fe] abundances of the most Rb-rich stars 
with [Rb/Fe] > 2. 

In regards to reaction rates, we tested different nuclear networks and different rates for 
the ^^Ne(Q;,n)^^Mg and ^^Ne(a,7)^^Mg reactions. Tests with a full s-process network of 320 
species shows that at solar metallicity the stars are not expected to produce many elements 
beyond the first s-process peak at Rb-Sr (except perhaps in the case when a ^^C pocket is 
included). In regards to the ^^Ne +a reactio ns, the NACRE rate s result in the highest Rb 



abundances. The newer ^^Ne +a rates from llliadis et al.l (120 lOf ). which have significantly 



lower uncertainties, produce [Rb/Fe] ratios that are lower by up to 0.3 dex, whereas the rates 
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from iKarakas et al.l (l2006bl ) produce final [Rb/Fe] that are up to 0.5 dex lower. Note that 



this is mostly because the Ne(a,7) Mg rate froni|I 



iadis et al. (12010) has been significantly 
reduced, while th e ^^Ne(a,n)^^Mg ra te is similar to iKarakas et al.l ( l2006bl ). Using the ^^Ne 
+a reactions from llliadis et al.l (|2010[ ) produces final maximum [Rb/Fe] ~ 0.9 dex and covers 
the lower end of the distribution of [Rb/Fe] in the observed OH/IR AGB stars. To match 
the composition of most of the Rb-rich stars requires us to use the faster NACRE rates. This 
may provide an important constraint on these notoriously difficult to estimate reaction rates. 
The main issue with this at the moment is the large uncertainties on the Rb abundances 
derived from observations. 

We finish with a discussion of the impact of the nucleosynthesis of the brightest AGB 
stars on the chemical evolution of galaxies and stellar systems. Delaying the superwind in 
intermediate-mass stars results in a far greater production of neutron-capture elements at 
solar metallicity and we would expect that this increased pr oduction to carry on through to 
the lowest metallicities. Our rnodel s and the observations by iGarcia-Hernandez et al.l (120061 ) 
and iGarcia- Hernandez et al.l (120091 ) show that the most massive C-0 core AGB stars of solar 
metallicity (up to 7Mq) and at the metallicities of the Large and Small Magellanic Clouds 
experience considerable dredge-up. This has some important consequences. If intermediate- 
mass stars at the metallicities of the Galactic globular clusters (GCs; —2.3 ^[Fe/H]< —0.7) 
also experience deep TDU, then intermediate-mass AGB stars would not be g ood candidates 
to ex plain the abundance anomalies observed in ev ery well studied GCs (e.g.. 
2004 ). This has been explored quantitatively by iFenner et al.l (120041 ) and 



Gratton et al. 



Karakas et al. 



(l2006al ) using yields from models with deep dredge-up, similar to those presented here. 
Intermediate-mass stars of up to 6 — 7Mq might instea d be candidates fo r producing the 
neutron-capture elements i n some GCs, includin g M4 ( Yong et al.l l2008bl Jal) and M22 (in 
particular see discussion in iRoederer et al.ll201ll ). Note that the neutron-capture elements 
observed in GCs do not show any correlation with the light-element abundance patterns, 
ruling out a relation between the two. However, models of stars with masses near the limit 
of C-0 core production (here the 8Mq model) or super-AGB stars with 0-Ne cores ma y 
not experience very efficient TDU (see the models of ISiessll2010l : IVentura &: D'Antonall201ll ). 
These objects experience very hot HBB and deper iding upon the r aodel, little pollution 



from He-intershell material (notably primary ^^C). ISiessI (120 lOl ) and IVentura fc D'Antona 
(120111 ) explored this idea using super-AGB models of the appropriate metallicity. Further, 
our results and theirs show that super-AGB stars can produce the high helium abundances 
of y ~ 0.4 needed to expl ain the multiple-p opulations observed in the color-magnitude 
diagrams of some GCs (e.g., iPiotto et al.ll2005l ). 



How can we test that intermediate-mass AGB stars of low metallicity did (or did not) 
produce the neutron-capture elements in GCs? Unfortunately it is not possible to verify 
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this production directly with observations today, ahhough there do exist bright s-process 



rich AGB stars in the SMC with metalhcities of [Fe/H] = —0.7 (JGarcia-Hernandez et aL 



2009f). whic h are at the upper end of the GC metalhcity distribution (e.g., 47 Tucanae, 



Harris! Il996[ ). However, the most metal-poor, massive AGB stars would have evolved away 
a long time ago. One hope is that future large-scale surveys of halo stars (e.g., the High 
Efficiency and Resolution Multi-Element Spectrograph for the AAT or SkyMapper) might 
reveal hints of such nucleosynthesis among the carbon and nitrogen-enhanced metal-poor 
star population. 
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Fig. 1. — The temporal evolution of the (a) growth of the H-exhausted cores, (b) the He-shell 
temperatures, and (c) the temperatures at the base of the convective envelope during the 
entire AGB. The black solid lines shows data for the 5Mq, Z = 0.02 model with a delayed 
superwind, and the red dashed lines show data for the 5Mq with IVassiliadis fc WoodI (Il993l ) 
mass loss, shifted to the left by t = 1 x 10^ years for clarity. (A color version of this figure 
is available on-line). 
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Fig. 2. — The temporal evolution of the (a) pulsation period over the first 40 TPs, (b) the 
growth of the H (black solid line) and He-exhausted (red dashed line) cores over the first 40 
TPs, (c) the He-shell temperature over the entire AGB, and (d) the temperature at the base 
of the convective envelope during the entire AGB for the OM©, Z = 0.02 model. 
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Fig. 3. — The evolution the stellar luminosity (in Lq) with TP number for the 9Mq, Z 
model. (A color version of this figure is available on-line). 
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Fig. 4. — The evolution of elements [X/Fe] lighter than Fe as a function of atomic number, 
Z, for the SM© model with a delayed superwind and nuclear network using 320 species. 
Included are the approximate locations (in proton number, Z) of some key elements. Each 
dot represents the surface composition after a thermal pulse. 
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Fig. 5. — The evolution of elements [X/Fe] heavier than Fe as a function of atomic number, 
Z, for the 5Mq model with a delayed superwind and nuclear network using 320 species. 
Included are the approximate locations of some key elements, noting that Pb is at Z = 82. 
Each dot represents the surface composition after a thermal pulse. 
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Fig. 6. — The evolution of elements [X/Fe] heavier than Fe as a function of at omic number, 



Z, for the QMq, Z = 0.0001 model with scaled solar initial abundances from iLugaro et al. 



(I2OI2I ). The nu clear network used for this calculation had 320 species and included the ^^Ne 
+ a rates from iKarakas et al.l (J2006bl ). Included are the approximate locations of some key 
elements, noting that Pb is at Z = 82. Each dot represents the surface composition after a 
thermal pulse. 
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Fig. 7. — The bolometric luminosity at the tip of the AGB versus the [Rb/Fe] abundance 
from the last computed TP (connected by the solid lines) and f rom the synthetic ev o lution 
calculations (connected by the dashed lines). Models using the IVassiliadis fc WoodI (119931) 
mass-loss prescription are connected by the lower solid and dashed lines, and models calcu- 
lated using a delayed superwind by the upper solid and dashed lines, respectively. Symbols 
and labels denote the initial stellar mass. The shaded region indicates the range of observed 
[Rb/Fe], noting that the maximum reaches +2.5 dex. The average observed [Rb/Fe] = 1.4 
is shown on the left along with the maximum uncertainty of ±0.8 dex. (A color version of 
this figure is available on-line). 
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Fig. 8. — The pulsation period from the stellar evolutionary sequence versus [Rb/Fe] from 
the post-processing nucleosynthesis calculations for the ^Mq and TMq models with a delayed 
superwind. (A color version of this figure is available on-line). 



Table 1: Characteristic parameters of the AGB models. 



M' 



M' 



M\ 



TPs NTPs 



ip 



An 



M, 



dredge 



-^bce 



T^max /TfPeaE j. ^-AGB-tip 

-^HcshcU ^^-'bol ^^-'bol 



(Me) 


(Mo) 


(Me) 


(days) 






(years) 




(Me) 


(MK) 


(MK) 






5.0 


1.787 


0.914 


500 


25 


17 


1.29(10^) 


0.95 


6.47(10-2 


) 57.3 


344 


-6.32 


-6.21 


5.0 


1.641 


0.756 


700 


72 


64 


1.18(10^) 


0.95 


2.22(10-1 


) 73.2 


354 


-6.51 


-6.29 


6.0 


1.673 


0.761 


500 


33 


25 


7.87(10^) 


0.93 


6.67(10-2 


) 76.5 


351 


-6.58 


-6.34 


e.o'^ 


1.583 


0.670 


500 


36 


25 


7.77(10^) 


0.95 


6.71(10-2 


) 87.2 


355 


-6.92 


-6.41 


6.0 


1.628 


0.691 


800 


158 


150 


4.77(103) 


0.92 


3.56(10-1 


) 84.1 


366 


-6.90 


-6.45 


7.0 


2.150 


1.181 


500 


42 


33 


4.67(10^) 


0.90 


5.32(10-2 


) 87.7 


364 


-6.83 


-6.48 


7.0 


2.115 


1.122 


800 


188 


178 


3.50(103) 


0.90 


2.70(10-1 


) 91.4 


378 


-7.04 


-6.60 


8.0 


2.621 


1.568 


500 


59 


49 


2.20(10-3) 


0.89 


3.95(10-2 


) 96.8 


378 


-7.12 


-6.64 


9.0 


3.208 


2.019 


500 


163 


155 


5.37(10^) 


0.77 


3.27(10-2 


) 111 


402 


-7.55 


-6.89 



Note: MK denotes the temperature in x lO^K, NTPs denotes the number of TPs with peak temperatures over 

300xlO^K, M^^^ denotes the peak AGB bolometric luminosity, and M^^^ ~ '^ denotes the AGB-tip bolometric 

luminosity, (a) Computed with a = 2.5. All other models have a = 1.86. 
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Table 2: Selected key elemental abundances in [X/Fe] (in dex) and isotopic ratios (by number) at the stellar surface at 
the end of the computed evolution for the nucleosynthesis models. Also included is the average mass fraction of helium 
in the stellar wind, <X(He)>. 
Mass P N^t 22Ne <X(He)> CjO ^^C/^'^G [N/Fe] [0/Fe] [Na/Fe] ^VFef^^MgT^^M^^^M^T^^^ 



5.0 


500 


166 


K06 


0.321 


1.204 


14.00 


0.505 


-0.037 


0.227 


0.037 


5.524 


4.638 


5.0 


500 


172 


NAC 


0.316 


1.192 


13.46 


0.509 


-0.034 


0.213 


0.057 


5.108 


4.213 


5.0 


700 


166 


K06 


0.332 


0.772 


11.14 


1.228 


-0.045 


1.229 


0.114 


2.404 


1.700 


5.0 


700 


172 


NAC 


0.331 


0.717 


11.49 


1.245 


-0.042 


0.308 


0.157 


2.013 


1.378 


5.0 


700 


172 


110 


0.331 


0.718 


11.43 


1.242 


-0.045 


0.306 


0.140 


2.467 


2.828 


5.0 


700 


320 


NAC 


0.334 


0.720 


11.13 


1.245 


-0.044 


0.311 


0.160 


1.981 


1.375 


6.0 


500 


166 


K06 


0.341 


0.678 


14.25 


0.869 


-0.055 


0.699 


0.053 


5.220 


4.103 


6.0 


500 


172 


NAC 


0.337 


0.647 


14.71 


0.879 


-0.052 


0.235 


0.062 


4.750 


3.677 


6.0 


500 


320 


NAC 


0.338 


0.633 


13.48 


0.884 


-0.053 


0.235 


0.062 


4.718 


3.682 


6.0 


800 


166 


K06 


0.384 


0.887 


8.440 


1.410 


-0.237 


1.466 


0.215 


1.679 


1.157 


6.0 


800 


172 


NAC 


0.371 


0.711 


9.046 


1.423 


-0.216 


0.343 


0.251 


1.067 


0.723 


6.0 


800 


172 


110 


0.371 


0.710 


8.950 


1.420 


-0.221 


0.358 


0.200 


1.263 


1.514 


7.0 


500 


166 


K06 


0.362 


0.582 


12.43 


0.897 


-0.100 


0.744 


0.050 


4.734 


3.996 


7.0 


500 


172 


NAC 


0.355 


0.545 


12.52 


0.908 


-0.096 


0.251 


0.044 


4.144 


3.510 


7.0 


800 


166 


K06 


0.396 


0.943 


7.588 


1.326 


-0.375 


0.985 


0.174 


0.714 


0.924 


7.0 


800 


172 


NAC 


0.383 


0.761 


8.005 


1.339 


-0.361 


0.097 


0.208 


0.360 


0.426 


7.0 


800 


172 


110 


0.383 


0.762 


7.989 


1.337 


-0.363 


0.121 


0.157 


0.386 


0.768 


8.0 


500 


166 


K06 


0.367 


0.440 


9.462 


0.903 


-0.117 


0.623 


0.030 


2.139 


3.186 


8.0 


500 


172 


NAC 


0.366 


0.395 


9.282 


0.921 


-0.127 


0.208 


0.026 


1.676 


2.620 


8.0 


500 


172 


110 


0.362 


0.399 


9.506 


0.914 


-0.114 


0.188 


0.022 


1.690 


3.474 


8.0 


500 


320 


NAC 


0.361 


0.391 


9.342 


0.914 


-0.108 


0.171 


0.027 


1.569 


2.545 



00 



9.0 
9.0 



500 
500 



166 

172 



K02 

NAC 



0.387 
0.375 



0.342 
0.300 



6.555 
6.370 



0.962 
0.975 



-0.131 
-0.136 



0.495 
0.085 



0.02 
0.028 



0.357 
0.197 



1.114 
0.631 



Note: Net denotes the network used for the calculation, 
and the initial ^^Mg/^^Mg and ^^Mg/^^Mg ratios are 7.8! 
the calculation where K06 is the rate from Karakas et al 



Ihadis et al.l fcOlOf ). 



The initial C/0 ratio = 0.55, the initial ^"^C/^'^C = 89.4, 
and 7. 17, respectively. Here ^^Ne denotes the rate used in 
fl2006bh . NAC is the NACRE rate, and ILIO the rate from 



Table 3: Selected key heavy element abundances [X/Fe] (in dex) at the stellar surface at the end of the computed 
evolution for the nucleosynthesis models. 



Mass 


P 


Net 


22Ne 


Co/Fc 


Cu/Fe 


Se/Fe 


Kr/Fe 


Rb/Fe 


[Sr/Fe 


Zr/Fe 


Ba/Fe 


5.0 


500 


166 


K06 


0.174 


0.092 


-0.02 


0.016 


0.030 


0.012 


-0.077 


- 


5.0 


500 


172 


NAC 


0.259 


0.222 


0.053 


0.064 


0.099 


0.035 


-0.063 


- 


5.0 


700 


166 


K06 


0.508 


0.446 


0.164 


0.191 


0.275 


0.128 


-0.016 


- 


5.0 


700 


172 


NAC 


0.581 


0.681 


0.555 


0.615 


0.746 


0.422 


0.228 


- 


5.0 


700 


172 


110 


0.558 


0.558 


0.338 


0.346 


0.457 


0.214 


0.058 


- 


5.0 


700 


320 


NAC 


0.584 


0.667 


0.574 


0.652 


0.763 


0.497 


0.332 


0.080 


6.0 


500 


166 


K06 


0.220 


0.146 


-0.002 


0.032 


0.052 


0.022 


-0.072 


- 


6.0 


500 


172 


NAC 


0.288 


0.293 


0.132 


0.139 


0.195 


0.078 


-0.038 


- 


6.0 


500 


320 


NAC 


0.287 


0.282 


0.133 


0.157 


0.198 


0.096 


0.045 


0.018 


6.0 


800 


166 


K06 


0.628 


0.671 


0.420 


0.492 


0.635 


0.387 


0.166 


- 


6.0 


800 


172 


NAC 


0.639 


0.870 


0.860 


1.026 


1.148 


0.807 


0.631 


- 


6.0 


800 


172 


110 


0.648 


0.775 


0.668 


0.764 


0.886 


0.563 


0.376 


- 


7.0 


500 


166 


K06 


0.216 


0.168 


0.016 


0.049 


0.074 


0.033 


-0.068 


- 


7.0 


500 


172 


NAC 


0.248 


0.296 


0.175 


0.201 


0.266 


0.116 


-0.011 


- 


7.0 


800 


166 


K06 


0.513 


0.581 


0.361 


0.436 


0.571 


0.340 


0.132 


- 


7.0 


800 


172 


NAC 


0.498 


0.751 


0.741 


0.904 


1.021 


0.693 


0.522 


- 


7.0 


800 


172 


110 


0.517 


0.678 


0.592 


0.696 


0.808 


0.504 


0.327 


- 


8.0 


500 


166 


K06 


0.197 


0.202 


0.054 


0.091 


0.132 


0.060 


-0.055 


- 


8.0 


500 


172 


NAC 


0.183 


0.300 


0.221 


0.284 


0.364 


0.172 


0.034 


- 


8.0 


500 


172 


110 


0.194 


0.268 


0.157 


0.194 


0.256 


0.112 


-0.012 


- 


8.0 


500 


320 


NAC 


0.195 


0.303 


0.242 


0.322 


0.397 


0.224 


0.133 


0.022 


9.0 


500 


166 


K06 


0.155 


0.246 


0.106 


0.164 


0.242 


0.116 


0.070 


- 


9.0 


500 


172 


NAC 


0.121 


0.317 


0.269 


0.363 


0.460 


0.225 


0.171 


- 



00 
to 
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Table 4: Details and results of synthetic evolution calculations. 



Mass 
(Me) 



P 

(days) 



TPs remaining 



AM 

(Mo) 



X(Rb) 



intcrshcU 



Final [Rb/Fe] 
(dex) 



5.0 


500 


3 


0.417 


5.5(10-^) 


0.38 


5.0 


700 


2 


0.412 


5.6(10-6) 


0.96 


6.0 


500 


3 


0.187 


1.1(10-6) 


0.41 


6.0 


800 


4 


0.193 


7.0(10-6) 


1.34 


7.0 


500 


6 


0.194 


2.0(10-6) 


0.54 


7.0 


800 


7 


0.150 


8.0(10-6) 


1.26 


8.0 


500 


15 


0.106 


3.1(10-6) 


0.80 


9.0 


500 


60 


0.031 


5.3(10-6) 


0.90 



Note: Intershell compositions are from the 172 species network which include the NACRE 

rate for the ^^Ne neutron source. 
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